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6 ABSTRACT 

CO 

We describe a non-parametric method to reconstruct gas density and 
temperature profiles of galaxy clusters from observations of X-ray surface 
brightness, emission-weighted temperature and the ^-parameter through the 
Sunyaev-Zel'dovich effect. This method is based on the inversion of the 
projected profiles under the spherically symmetric approximation, which is 
independent of the equation of state of cluster gas and does not assume that the 
gas is in hydrostatic equilibrium. In particular we examine the reliability of the 
reconstruction method assuming a few theoretical models for cluster gas and 
assigning the statistical errors expected for future X-ray and radio observations. 
We also discuss briefly the effect of non-sphericity on the basis of simulated 
clusters. 
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1. Introduction 

The importance of cluster abundances and baryon fractions is well-recognized and they 
are now regarded as one of the standard established statistics to constrain the cosmological 
models (e.g., Henry & Arnaud 1991; White, Efstathiou, & Frenk 1993; White et al. 1993; 
Viana & Liddle 1996; Eke, Cole, & Frenk 1996; Kitayama & Suto 1997; Kitayama, Sasaki 
& Suto 1998; Shimasaku 1998). On observational sides, the accuracy and reliability of 
such statistics will be considerably improved with several upcoming projects. In particular 
the angular resolutions of the AXAF (Advanced X-ray Astrophysics Facility) and XMM 
(X-Ray Multi-Mirror Mission) are 0.5" in 0.1-10 keV and 15" in 0.1-15 keV, respectively, 
and therefore they are expected to provide high-resolution X-ray surface brightness and 
temperature maps for many clusters with unprecedented angular and energy resolutions. 
The existing interferometric facilities in radio bands including BIMA (Berkeley Illinois 
Maryland Association) and OVRO (Owens Valley Radio Observatory) (e.g., Cooray et 
al. 1998) also produce two-dimensional maps of clusters via the Sunyaev-Zel'dovich (SZ) 
effect with 10" ~ 30" resolution and signal-to-noise ratio more than 20. Furthermore, the 
Japanese Large Millimeter and Submillimeter Array (LMSA) project, for instance, will start 
an extensive survey of clusters in millimeter and submillimeter bands (see also Komatsu 
et al. 1998). Such detailed information of intracluster gas in multi-bands will significantly 
improve the role of clusters as cosmological probes. 

On theoretical sides, recent N-body/hydrodynamical simulations (e.g., Couchman 
et al. 1995; Bryan & Norman 1998; Eke & Frenk 1998; Yoshikawa, Itoh, & Suto 1998) 
start to reveal several important keys to understanding the intracluster gas on a physical 
basis. In particular, Navarro, Frenk & White (1997; NFW hereafter) suggest that a dark 
matter halo of cluster is described by a universal density profile, or at least by a class of 
well-specified functional forms (see Fukushige & Makino 1997; Moore et al. 1998). Their 
results imply that with appropriate equation of state, the gas density and temperature 
profiles of clusters can be predicted given a cosmological model, at least in a spherically 
symmetric approximation (Makino, Sasaki & Suto 1998; Suto, Sasaki, & Makino 1998). In 
other words, one does not have to adopt too simplified empirical models like the isothermal 
/9-model in approximating the cluster gas. In fact we are now in a position to construct a 
more physical model combining the new observational data and the theoretical predictions. 
By reconstructing the cluster profiles in a model-independent manner from the observations 
and comparing them with the theoretical models, one can directly test the universal dark 
halo conjecture, and if confirmed, probe the cosmological parameters in an independent 
and complementary manner with more conventional approaches like cluster abundances, 
cosmic microwave background anisotropies, and galaxy clustering statistics (e.g., Tegmark 
et al. 1998). Therefore a reliable and model-independent reconstruction of radial profiles 
of cluster gas should also contribute significantly to the conventional analyses of clusters 
including the estimates of the Hubble constant (Silk & White 1978; Inagaki, Suginohara, 
& Suto 1995; Kobayashi, Sasaki & Suto 1996) and peculiar velocity (e.g., Yoshikawa, Itoh, 
& Suto 1998), baryon fraction and gravitational lensing (e.g., Wu & Fang 1997). See 
Birkinshaw (1998) for an excellent review on the SZ effect and its cosmological implications. 
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In this paper, we apply an idea proposed earlier by Silk & White (1978), and examine 
in detail a non-parametric reconstruction method of the radial profiles of clusters, which 
does not depend on any assumption such as hydrostatic equilibrium and equation of state 
of cluster gas except for the spherical symmetry of clusters. We also propose a new method 
to reconstruct the temperature radial profile of clusters only from X-ray observations. This 
method is also applicable to other elliptic models as well, provided that the ellipticity 
is assumed a priori (c.f., Fabricant, Rybicki & Gorenstein 1984; Hughes & Birkinshaw 
1998). In this context, we note that Zaroubi et al. (1998) recently proposed a different 
reconstruction method of 3-dimensional profiles of clusters using the Fourier Slice Theorem, 
and our current scientific motivation is quite similar with theirs. While their method, 
employing Fourier transform of the projected images and inverse transformation, is fairly 
complicated and it is not clear yet to what extent it is practical, our procedure is rather 
straightforward and useful. 



2. Abel's integral solution for radial profiles of gas density and temperature in 

clusters of galaxies 

Most observable quantities for a cluster are often written as an integration over the 
line-of-sight: 

/oo rco vdr 

g{r)dl = 2 g{r) (1) 
-oo Jd A e J r 2 _ (p A 02 

where the cluster is assumed to be spherically symmetric, 9 and r denote the (projected) 
angular separation and (3D) spatial radius from the cluster center, and 6a is the angular 
diameter distance to the cluster. Using Abel's integral, equation (|l|) can be inverted to give 

g{r) = j_ r ,dA m d i (2) 

ird A Joc d6 ^e 2 -r 2 /d\ 

The above Abel's integral solution can be readily applied to the bolometric X-ray 
surface brightness: 



/oo i />oo 

a x (T e )n 2 e (Je*d 2 A (z) + P) dl = 2A X / a x (T e )n 2 e (r) - 
-oo Jdjsfi /. 



rdr 



r 2 - d 2 A 6 2 



(3) 



and y-parameter, or more generally, the SZ flux at a given frequency band: 

rdr 



rOO VflT 

S y {9)=2A y T e {r)n e {r) (4) 
Jd A e J r 2 _ ^2 



In the above expressions, T e (r) and n e (r) are the temperature and number density of the 
electron gas of the cluster, and ax(T e ) is the X-ray emissivity. While the coefficients A x 
and A y can be explicitly computed once the observational bands are specified, their specific 
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values are not important in the following analysis. Silk & White (1978) were the first to 
propose to apply the Abel inversion to the X-ray and SZ profiles, but their discussion is 
limited to the estimation of the central values of gas density and temperature aiming at 
the determination of (Ia{z). The purpose of the present paper is to examine the feasibility 
of the reconstruction method quantitatively, properly taking account of the current and 
upcoming data quality 

If one considers only the thermal bremsstrahlung for the X-ray emissivity, which is a 
good approximation for clusters with T e £ 3keV, then ax(T e ) = T e 1//2 (r), and equations (|3|) 
and (F§) are inverted to yield 



[n e (r)f 



A,. 



dS x {d) 



de 



de 



^e 2 - r 2 /d 



it Aid a. 



•/d A dS y (6) 



de 



(5) 



de 



[T e (r)] 2 /a(T c ) = [T e (r)] 3/2 



.4, 



e 2 - r 2 jd\ 

r/dA dS y (e) 



de 



de 



sjff 2 - r 2 /d 2 



nA 2 y d A 



r ' d * dS x (9) 



de 



(6) 



de 



e 2 -r 2 /d\ 



In reality, the observed X-ray flux is band-limited and one should also take into account 
other emission processes than the thermal bremsstrahlung. While the present methodology 
works for those cases, the result cannot be expressed in a simple form as the above equation, 
and we focus on the case of otx(T e ) = T*' 2 (r) as an explicit example. 
If the emission-weighted temperature profile projected on the sky: 



T e a x (T e )n 2 e ^d 2 A e 2 + P) dl f°° T c a x (T c )n 2 c ^d 2 A e 2 + I 2 

J — CO J — CO 



dl 



T X (9) 

/ ax(T e )nl^d\e 2 + l 2 ) dl 

J — oo 

is measured, however, one can similarly derive 
T e (r) « X (T C ) [n e {r)f 



s x (e)/A x 



(7) 



•IdA d[S x (6)T; 
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Combining with Abel's integral solution of equation 
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(9) 
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one can obtain the temperature profile T e (r) by 

r/dA d[S x {9)Tx{6)] d6 



de 



e 2 - r 2 /d\ 



Te(r) = T^Tds:, Te • (10) 



de Je 2 - (r/d A ) 2 



Equation |X0| is more suitable for reconstructing the temperature profile than equation (y) 
in that all the necessary information can be obtained from the X-ray observation alone, 
and that it is independent of the specific form of the X-ray emissivity ax(T e ). The use of 
equation (^) requires sufficiently good spatial and energy resolutions, both of which are 
feasible only with future X-ray satellites including AXAF and XMM. 



All the above procedures (eqs.||, @ and [K|) involve the integral of the form: 

r / d A df{6) de 

de Jm-{r/d A ) 2 ' 



In what follows we numerically evaluate the integral as 

Jmax -P -f 

— +1 ('\D 

where the observed quantity fi (i = l,i max ) is given at the discrete angular radius 9i, and 
e i+ i = [6i + e i+ i)/2. i max is the index for the outermost bin and i m i n indicates the bin 
which corresponds to r/dA- In practice we adopt i max = 30 for definiteness. In fact one 
could improve the above evaluation by pre- smoothing the dataset to reduce the discrete 
noise component. While such a procedure is inevitably model-dependent and we did not 
attempt one in what follows, this would be practically useful in dealing with the actual 
observational data and improve the results below. 



3. Statistical and systematic errors in the reconstruction 
3.1. models of gas density and temperature profiles 

Strictly speaking one needs to observe S x (e) and S y (9) for r/d^ < 6 < oo in order to 
determine n e (r) and T e (r). In practice, however, the observable fluxes are limited to the 
finite extent of clusters, and also contaminated by the angular resolution of the telescope 
and the detector noise especially at outer regions. Moreover realistic clusters are not 
spherically symmetric to some degree. In this section, we employ several cluster models 
with different density and temperature profiles, and examine how the above reconstruction 
method works in taking account of such observational limitations. More specifically, we 
examine the following four models: 
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Model A: this is a conventional isothermal /?-model described by 

= [1 + (r / e r ° c)2] 3^/ 2 ^ W = T eO- (13) 

We adopt a typical value for clusters, (3 = 2/3, for definiteness. The values of the other 
parameters, n e0 , T e0 , and r c , are irrelevant for our current purpose since all the results can 
be rescaled appropriately. 

Model B: consider a family of the dark matter halo profiles given by 

Pom(x) = x , { ^ xU)X , (14) 

where x = r/r s is the dimensionless radius in units of the characteristic scale r s , and 5 DM is 
the amplitude of the profile. NFW claimed that a halo profile with (//, z/, A) = (1, 2, 1) is the 
universal shape fairly independent of the cosmological initial conditions (see also Fukushige 
& Makino 1997 and Moore et al. 1998). 

Suto, Sasaki & Makino (1998) found an analytical solution for the gas density 
and temperature profiles when the gas obeys the polytropic equation of state 
T e (r) = T e o[n e (r)/n e o] 1//n and its self-gravity can be neglected compared with the 
dark halo. Their solution is written as 

n e (r) = n c0 [l - B p f(x)} n , T e (r) = T c0 [l - B p f(x)}, (15) 

where 

rx U^~^ i f x U^~^ 

fix) = / — -^du / — -r-rf-u, (16) 

Jy ' Jo (1 +u u ) x xJo (1 + u») x ' v ; 

B = AnG ^9 m p 5 ™ r l (17) 
p ~ n + 1 kT e0 ' 1 ' 

G is the gravitational constant, k is the Boltzmann constant, fi g is the mean molecular 
weight of the gas, and m p is the proton mass. 

Therefore given the dark matter halo profile, the gas density and temperature profiles 
are specified by the additional two parameters n and B p . We consider the following two 
models (model B and C) which have the analytic solutions for f(x). In the NFW model 
(/i — 1, v — 2, and A = 1), 

/w = 1 _Ml±f). (18) 

x 

We adopt n = 12, B p — 1.0 for model B so that the resulting S x (9) has a fairly similar 
profile to that of model A (corresponding to (3 = 2/3). 
Model C: halo model with /i = 3/2, v — 1, and A = 3/2 has 



f( x ) = 2 J - - In (v^ + VT+x) ■ (19) 
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We adopt n = 11, -B p = 0.5, again so that the resulting S , a; (6') resembles the profile of model 
A. 

Model D: we select a simulated cluster (cluster A at z = in Yoshikawa, Itoh, & Suto 
1998) in order to examine a possible effect of the non-sphericity. The cluster has lower 
temperature in the central region unlike the above three models; in this sense also, it can be 
regarded as a good example to test the robustness of the present reconstruction method. In 
practice, we find that the circularly averaged profiles of its X-ray and SZ surface brightness 
of the cluster are well approximated by 

Sx{6) = [l + {d A 9/r x y]Wx-W S ^ 9) = [l + (d^/r[)2]3ft/2-i/2' ( 2 °) 

with p x = 0.69, r x = 0.13Mpc, (3 y = 0.86, and r y = 0.18Mpc. Therefore we adopt the 
above spherical fits (with appropriate errors; see the next subsection) as input data, and 
examine the extent to which the reconstructed density and temperature profiles agree with 
the spherically averaged ones directly computed from the simulated data. 



3.2. error assignment and reconstructed profiles 

Figure [1] shows the input projected profiles of bolometric X-ray flux, y-parameter and 
X-ray emission- weighted temperature (eqs. ||, Q, and 0) for models A to D described 
in §3.1. In order to estimate statistical and possible systematic errors to our method due 
to observational uncertainties and deviation from spherical symmetry, we perform the 
reconstruction after adding the following Gaussian distributed errors. More specifically we 
examine models A, B, and C to estimate errors from observational uncertainties, and model 
D from the deviation from spherical symmetry. 

The relevance of the error assignment should sensitively depend on both the observing 
facilities and the specific target clusters one has in mind. For models A to C, we consider 
two cases (I and II) for error assignment; for S x the assumed errors are fairly realistic even 
with the current X-ray satellites (e.g., Briel & Henry 1996) at least for relatively near 
and bright clusters, but for S y and T x the errors should be feasible only with the future 
interferometer facilities including Japanese LMSA (Large Millimeter and Submillimeter 
Array; Kawabe, private communication), and future X-ray satellites with high sensitivity 
such as XMM with sufficiently large integration time. For a given error of S x , one may 
roughly scale the resulting error-bars in Figure ^ for the different level of the errors in S y 
and Tx- With this in mind, we consider two specific cases; for the case I, we assign the la 
errors as 

AS X {9) = 10-^(0), AS y (6) = 10- 2 ^(0), AT X (6) = 10- 2 T X (9) (21) 
at each 9, and for the case II they are 

AS X (9) = 10~ 3 S X (0), AS y (9) = 2 x 10^(0), AT X (9) = 5 x 10~ 2 T X (#). (22) 
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Error bars in Figure [I] indicate those la dispersion of Gaussian errors in the case of I and 
II (all the curves for case I are artificially shifted upward so that the two cases are easily 
distinguished). For model D, on the other hand, we assign the la dispersion of the fluxes 
around the circularly averaged values directly computed from simulation data profile by 
Yoshikawa, Itoh, & Suto (1998), and do not attempt to include the observational errors. 
We sample the data using 30 radial bins in logarithmically equal intervals. 

According to the above procedure, we construct 200 realizations for each model with 
different random number for the error assignment, and then perform the reconstruction. 
Figure shows the reconstructed profiles of gas density (top panels) , the ratio between 
them (upper middle panels), temperature (lower middle and bottom panels) for each 
model. The lower middle panels show the temperature profiles reconstructed from X-ray 
surface brightness and the SZ flux (eq. ||) and the bottom panels show that reconstructed 
from X-ray surface brightness and emission- weighted temperature (eq. |10|). Solid lines 
correspond to the true profiles for models A to C (again the results for case I are artificially 
shifted upward for an illustrative purpose), and the spherically averaged profiles for model 
D. 

The reconstructed profiles are plotted with their la error-bars computed from the 
200 realizations. Since the main contribution for the projected fluxes comes from the 
radius around r c , the estimates become less reliable either for r < r c or r > r c . Also 
the reconstruction works better for n e (r) than for T e (r) since the X-ray emissivity is 
more sensitive to the former (cx n 2 e T^ 2 for the thermal bremsstrahlung) . The upper and 
lower solid lines in model D represent the ±la deviation from the spherically averaged 
profile, which illustrate the degree of non-sphericity of the simulated cluster. The degree of 
asphericity in our simulated cluster is comparable to, or even larger than, the error due to 
the reconstruction procedure. Since our simulated cluster seems to be typical in the light 
of the distribution of non-sphericity in the observed cluster sample (Mohr et al. 1995), the 
non-sphericity of that degree does not seriously degrade the current methodology although 
the inclusion of non-spherical effect is definitely an important next step (Zaroubi et al. 
1998). 

4. Discussion and conclusions 

We have presented a method to reconstruct the radial profiles of gas density and 
temperature in clusters of galaxies. Since most existing techniques start with some sort of 
empirical cluster gas profile and attempt to find the best-fit in that modeling, the result is 
fairly model-dependent. Our current method is in marked contrast to the previous ones in 
that the radial profiles can be reconstructed in a non-parametric manner without assuming 
the equation of state or hydrostatic equilibrium of cluster gas if high-resolution projected 
profiles of X-ray and SZ observations are given. 

We apply this method to three analytic spherical models (A to C) with gas density 
and temperature profiles and one non-spherical simulated cluster. For models A to C, the 
fractional uncertainties of reconstructed gas density profiles are less than 1% in our case I 
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while they amount to 20-30% at central regions and 2-3% at outer regions in our case II. 
On the other hand, the estimate of temperature is sensitive to the accuracy of the SZ flux 
or X-ray emission weighted temperature as in Figure 0. With the degree of asphericity in 
our simulated cluster, the resulting systematic error is smaller than our adopted statistical 
errors. In fact while our error assignment for X-ray temperature and the SZ fluxes may be 
rather optimistic, one can rescale the resulting error bars according to the real observation 
if necessary. 

Although we have presented the analysis on the basis of the spherical symmetry, the 
same methodology is in principle applicable to other non-spherical systems if one accepts 
some additional assumption on the density distribution along the line of sight. In this 
context, one of the most practically important ways is to assume a bilateral symmetry (e.g., 
Fabricant, Rybicki & Gorenstein 1984; Hughes & Birkinshaw 1998). In this case, all the 
results presented in this paper is readily applicable simply by an appropriate choice of the 
coordinate transformation. 

Considering the future extensive observation of clusters of galaxies in X-ray and 
radio bands with upcoming facilities such as ABRIXAS, AXAF, XMM, and LMSA, 
high-resolution images of clusters of galaxies in X-ray and radio bands will be available 
with unprecedented high signal-to-noise ratio and the reconstruction of profiles of clusters 
of galaxies with sufficient accuracy should be feasible. More quantitative discussion of 
the error and the resulting accuracy of the reconstruction should definitely depend on the 
specific target cluster, and we are not able to proceed further at this moment. Nevertheless 
we believe that the current methodology becomes useful in near future. Furthermore 
with such a reliable reconstructed radial profiles of clusters of galaxies, even if feasible 
only for relatively near and bright ones, one can revisit various cosmological issues, in 
which simplified and/or empirical models have been adopted, including the estimates of 
Hubble constant (Kobayashi, Sasaki & Suto 1996), peculiar velocity of clusters of galaxies 
(Yoshikawa, Itoh, & Suto 1998), the L — T relation which is inconsistent with observations 
if derived from a simple scaling argument, and physics cooling flow. 

We thank Pat Henry for information of the current observational errors of the X-ray 
fluxes, and Ryohei Kawabe for useful discussion on the current and future observing 
facilities in radio, millimeter and submillimeter bands and their expected performance. We 
also thank an anonymous referee for several comments which helped improve the earlier 
version of the paper. Numerical computations were carried out on VPP300/16R and 
VX/4R at AD AC (the Astronomical Data Analysis Center) of the National Astronomical 
Observatory, Japan, as well as at RESCEU (Research Center for the Early Universe, 
University of Tokyo) and KEK (National Laboratory for High Energy Physics, Japan). 
This research was supported in part by the Grants-in-Aid by the Ministry of Education, 
Science, Sports and Culture of Japan (07CE2002) to RESCEU, and by the Supercomputer 
Project (No. 97-22) of High Energy Accelerator Research Organization (KEK). 
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[Input projected profiles] 




e/e c e/e s e/e s e/e c 



Fig. 1. — Input profiles of scaled X-ray surface brightness (upper panels), the SZ flux 
(middle panels) and emission- weighted temperature (lower panels) for models A to D. For 
models A, B, and C, the smaller and larger error bars correspond to our cases I and II of 
error assignment, respectively (eqs. [[ZT|] and [22]). The curves corresponding to case I are 
artificially shifted upward for an illustrative purpose. For model D, the error bars indicate the 
la deviation from spherical symmetry obtained from the numerical simulation (Yoshikawa 
et al. 1998). The angular core radius 9 C in model D is computed from the core radius r c 
of the /3-model fit to the spherically averaged density profile of the simulated cluster. Solid 
lines indicate the true profile before assigning the error. 
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[Reconstructed radial profiles] 




10" 1 1 10" 2 10-' 1 10" 2 10-' 10" 1 1 10' 

r/r c r/r 3 r/r s r/r c 



Fig. 2. — Reconstructed gas density and temperature profiles for models A to D. Top panels: 
reconstructed n e (r) in units of the central value; Upper middle panels: reconstructed n e (r) in 
units of the true density n e)true (r); Lower middle panels: T e (r) reconstructed with X-ray and 
SZ fluxes in units of the central value; Bottom panels: T e (r) reconstructed with X-ray flux 
and emission-weighted projected temperature in units of the central value. For models A to 
C, solid lines indicate the true model profile, and the error bars correspond to those shown in 
Figure 1. The curves corresponding to case I are artificially shifted upward for an illustrative 
purpose. For model D, the middle solid lines in each panel represent the spherically averaged 
profiles, and the upper and lower solid lines indicate the ±ler dispersion of the simulated 
cluster around the spherically averaged profile. 



